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I. INTRODUCTION 


The research goals of this project are to provide physical insight 
into the nature of the extended gaseous atmospheres of both outer planet 
satellites and comets. For the outer planets, research efforts are fo- 
cused upon understanding the large circumplanetary atomic hydrogen dis- 
tribution in the Saturn system, both in terms of its neutral sources 
(which include the satellite Titan) and its role as a plasma source for 
the planetary magnetosphere. For comets, the emphasis is to understand 
the basic chemical composition of the nucleus and its interaction with the 
sun and solar wind through study of the composition and spatial structure 
of its extended cometary hydrogen, oxygen and carbon atmospheres. 

To understand the circumplanetary atomic hydrogen in the Saturn sys- 
tem, the overall strategy has been to model the spatial structure of the 
hydrogen torus of Titan (the expected dominant source of H atoms) and 
compare it with the Voyager UVS data for the hydrogen Lyman-a emissions. 

To achieve this objective, a collaborative effort has been established 
with D.E. Shemansky to provide optimally prepared Voyager UVS Lyman-a data 
for this comparison. The Titan hydrogen torus model at AER has also been 
improved to include the spatial lifetime of H atoms in the planetary mag- 
netosphere by incorporating the best available Voyager PLS electron and 
ion data. This model-data comparison will not only allow the Titan source 
to be quantified but will also provide a means of identifying and as- 
sessing the relative importance of other possible non-Titan hydrogen 
sources (i..e., the icy satellites, the planetary rings, and the planetary 
atmosphere) . 

The basic chemical composition of the comet nucleus, which is too 
small to be seen and is furthermore obscured from view by the gas and dust 
coma, will be investigated by studying the observed nature of the 
very-extended atmospheres of the comet. For this strategy (which is wide- 
ly adopted) to be successful, however, it is vital to have a model for the 
cometary atmosphere that accurately contains all the relevant physical 
interactions that occur in the sun/solar wind environment for the 
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so-called parent molecules that are ejected from the nucleus. Recognizing 
this requirement, a very general particle trajectory model (PTM) has been 
developed at AER. In this project, this model is used to analyze the den- 
sity and UV emissions of extended cometary hydrogen, oxygen and carbon at- 
mospheres. For this purpose, a collaborative effort with A.l.F. Stewart 
has been established to provide cometary H, C, and 0 emission data 
obtained by the ultravoilet spectrometer of the Pioneer Venus Orbiter for 
Comets P/Encke, P/Giacobini-Zinner , and P/Halley. 

The discussion of third year progress and achievements is emphasized 
in this final report, with a brief overview summary of first and second 
year efforts provided for completeness. A more thorough discussion of the 
accomplishments in the first and second year of this project is contained 
in the 1985 and 1986 annual reports. 
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II. HYDROGEN DISTRIBUTION IN THE SATURN SYSTEM 


1 . Overview 

The initial three-year plan for this project in the Saturn system is 
summarized in Table 1. Although the basic strategy has remained the same, 
cometary emphasis because of the return of P/Halley and a discovery of 
significant amounts of new Voyager data for the circumplanetary distribu- 
tion of hydrogen at Saturn have altered the implementation of this plan as 
summarized below. 

In the first year, scheduled objectives were achieved for refining 
the lifetime description of hydrogen atoms in the magnetosphere and for 
evaluating and partially processing Voyager UVS data for the Lyman-a emis- 
sion. A new lifetime process for hydrogen involving collisional loss with 
H atoms of the interstellar medium was identified. Preliminary model cal- 
culations were performed to determine the best procedure for including the 
long lifetimes of hydrogen atoms. 

In the second project year, efforts to have been expended on Saturn 
system research were largely postponed until the third year. This was be- 
cause of the major emphasis placed on cometary activities with the return 
of Comet P/Halley and because of difficulties in acquiring adequate 
self-consistent plasma information from the Voyager PLS instrument from 
which to calculate the hydrogen atom lifetime in the planetary 
magnetosphere . 

In the third year, the program outlined in Table 1 was resumed. The 
primary development, however, has been centered about the discovery by 
D.E. Shemansky in our collaborative effort of a significant amount of new 
Lyman-a data for hydrogen in the old Voyager UVS data tapes. This dis- 
covery, which is discussed below, has established the basis for a whole 
new set of studies and suggests that both Titan's hydrogen torus and 
Saturn's hydrogen corona play a major role in the circumplanetary dis- 
tribution of hydrogen. The new studies have provided the basis for re- 
search that is to be undertaken in a new NASA- sponsored project that 
logically continues the efforts initiated here. 
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SATURN SYSTEM: THREE-YEAR PLAN FOR MODELING ANALYSIS 
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2. Hydrogen Lifetime in the Saturn System 

The four relevant loss processes for atomic hydrogen in the Saturnian 
system are summarized in Table 2. The first two processes require a spa- 
tial description of the plasma properties in the planetary magnetosphere 
and also in the solar wind beyond the magnetosphere. The plasma descrip- 
tion for the solar wind is readily available. The plasma description for 
the magnetosphere is more difficult to specify because it must be deter- 
mined from limited analysis of data acquired by the Plasma Science (PLS) 
experiment on the Voyager 1 and Voyager 2 spacecrafts during their en- 
counters with the planet (12 November 1980 and 26 August 1981, respective- 
ly). The third loss process of Table 2, photoionization, is easily evalu- 
ated. The fourth process, elastic collision with atomic hydrogen in the 
interstellar medium, was identified last year by this project. Evaluation 
of this loss rate requires the determination of low-velocity 
cross-sectional information currently unavailable. This cross-sectional 
information for hydrogen atoms is being calculated by Shemansky as part of 
his collaborative involvement with this ongoing project. 

To construct a plasma description for the planetary magnetosphere, 

PLS data acquired by the Voyager spacecrafts have been utilized. The PLS 
data is divided into the reduction and analysis of electron data (being 
performed at the Goddard Space Flight Center, primarily by E.C. Sittler) 
and ion data (being performed at the Massachusetts Institute of 
Technology, primarily by J.D. Richardson). As noted in the second year 
annual report, significant discrepancies exist between the Voyager PLS 
electron and ion data for the magnetosphere. In addition, large and like- 
ly unrealistic differences exist between the amount of plasma inferred on 
L- shells in the inbound and outbound path of the Voyager spacecraft. To 
properly address these issues, E.C. Sittler has agreed to work with J.D. 
Richardson and publish a paper resolving these matters. Prior to this 
publication, their results will be available to the continuation of this 
project in which modeling analysis will be undertaken for the complete set 
of Lyman-a emission data for the Saturn system obtained by the UVS 
instrument of the Voyager 1 and 2 spacecrafts during their 1980 and 1981 
encounters with the planet. 
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Table 2 


LOSS PROCESSES FOR ATOMIC HYDROGEN IN THE SATURNIAN SYSTEM 


1, H + 


e (cold) 
e (hot) 


► H + + 2e 




s > 

H + 


s > 

H 

2. 

H + 

0 + 

► H + + ■ 

0 • 



N + 


N 



V J 


^ 4 


3 , H + hv 


► H + + 


4. H + H 


► H + H 


ISM TORUS 

(~20 km sec -1 ) (slow) 


ISM TORUS 

(fast) (fast) 
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3. Voyager UVS Data for Atomic Hydrogen 

Voyager measurements published to date for the hydrogen Lyman -a 
distribution in the Saturn system are for radial -scan data obtained in the 
Voyager 1 pre - encounter period (Broadfoot et al. 1981) and vertical-scan 
data obtained in the Voyager 2 pre-encounter period (Sandel et al. 

1982). As part of a collaborative effort with this project, D.E. 

Shemansky has prepared additional Voyager data for our modeling analy- 
sis. In the second project year, significantly improved radial scan data 
were identified and reduced by Shemansky for the Voyager 2 pre -encounter 
period. Additional data of lesser quality were also examined in hopes of 
constructing a complete image of the Lyman-o distribution. In this 
period, a search of the magnetic tapes containing all pertinent UVS data 
of Voyager 1 and Voyager 2 was undertaken. The search has uncovered a 
significant amount of new data that had not previously been reduced. In 
fact , this new data set is approximately 2 . 5 times larger than the old 
data set that was already reduced and published by the UVS team. A sum- 
mary of the older published data and newly reduced data is given in 
Table 3. 

There are a number of interesting aspects of the new data in Table 3 
that deserve discussion. The most spectacular is that from the Voyager 1 
post -encounter data from which an image of the H Lyman-a emission has been 
constructed at a viewing angle of -25® to the orbit plane of Titan. A 
preliminary version of this image is shown in Figure 1. The image indi- 
cates that the brightness distribution is not cylindrically symmetric as 
expected for a Titan torus alone (Smyth 1981; Hilton and Hunten 1988). 

The brightness distribution between Saturn and Titan's orbit has a minimum 
at pre -dawn and a maximum in the vicinity of the dusk terminator line sug- 
gesting an asymmetric dayside hydrogen source for a large electroglow- 
driven planetary corona. In the less corona -dominated portion of the cir- 
cumplanetary hydrogen, the presence of the Titan hydrogen torus can be 
distinctly identified. The image data are discussed in detail by 
Shemansky and Hall (1988). In addition to these data, the newly reduced 
Voyager 2 pre -encounter data provide the best one -dimensional data in the 
orbit plane. Combining the earlier published data and new data in the 
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Table 3 

Voyager Spacecraft Data Available for Research on the Atomic 
Hydrogen Distribution in the Saturn System 


Vovager 1 pre- encounter data 

Published Data: 14 days of one -dimensional scan data nearly 

parallel to the satellite orbit plane have been reported in the 
VI 30 day science report (Science 212 . 201-211, 1981). These 
data were obtained at a range of - 2 x 10 7 km from Jupiter. 

New Data: Additional data of similar nature is available but 

the quality is uncertain. The reduction process should be 
repeated to include the additional data. 

Voyager 1 post-encounter data 

Published Data: none 

New Data: Observations were obtained in mosaic scans across the 

system. All the data have been reduced for the first time in 
the past few months of 1987. The Voyager 1 post encounter 
trajectory is out of the solar system plane, and the data are 
obtained at an angle of -25° to the orbit plane of Titan. An 
image of the system in H Lyman-a emission has been constructed 
but is not yet published. The data in this image were obtained 
over the period 1980 DOY 324-343 (6-25 days post encounter) at 
spacecraft planet ranges of 0.83 x 10 7 - 0.33 x 10 s km. The 
integration time of the data set is 126 hours. Titan moves 
-1.25 orbits during the 20 day observing period. 

Voyager 2 pre -encounter data 

Published Data: Observations of one -dimensional scan data ap- 

proximately normal to the satellite orbit plane were published 
in the V2 30 day science report (Science 215, 548-553, 1982). 

New Data: All usable observations for one -dimensional scan data 

approximately parallel to the satellite orbit plan have been re- 
duced for the first time. The data were obtained in 1981 DOY 
180-186, at ranges of 0.54 x 10* - 0.49 x 10* km. The in- 
tegration time for this data set is 88 hours. These observa- 
tions provide the best one -dimensional data in the orbit plane 
of Titan that were obtained by the Voyager spacecrafts. 

Voyager 2 post -encounter data 

Post -encounter sequences were not obtained because of the lockup 
of the scan platform drive. 
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Figure 1 

Lyman -a Image for Hydrogen in the Saturn System 
(see page 15 for legend) 
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Figure 1. Density plot of the Saturn system in H Ly a emission obtained 
from Voyager 1 UVS mosaic scans of the system in the period, 

1980 DOY 324-343 (6-25 days post encounter). Spacecraft -planet 
range 0.83 x 10^ - 3.3 x 10^ km. The total integration time for 
the observations was 126 h. The pixel size is 1 Rj x 1 Rj 
perpendicular to the spacecraft-planet line, 1 Rj x 2.4 Rj 
projected onto the satellite orbital plane. The zenith angle of 
the spacecraft is 65° . The data is shown after subtraction of 
the LISM background calculated with the aid of a model, and weak 
stellar sources which appear in three locations approximately 
along the right side of the major axis of the Titan orbital 
ellipse drawn on the figure. The bright spot at the right ansa 
is due to imperfect correction for one of the stellar sources. 
The location of the sun-planet line is indicated in the upper 
right sector. The planet is drawn to size on the figure with 
the superposed mean size and orientation of the Voyager UVS 
field of view. A small sunlit crescent affects the signal on the 
right side of the planet. The data has been smoothed in a 3 x 3 
pixel area, weighted by integration time, in a 1,2,1; 2,4,2; 
1,2,1 pattern. The brightness scale in Rayleighs is given on 
the figure. The depression in signal in a rectangular area 
oriented parallel to the length of the slit at the position of 
the planet is caused by occultation of the LISM signal by the 
planet. The darkside of the planet would otherwise be of the 
same order of brightness as the coronal hydrogen. Bright areas 
toward the periphery and in some locations in the vicinity of 
Titan's orbit on the left side lose their significance because 
of low signal to noise ratios. 
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analysis of the H distribution with an appropriately developed model for 
the planetary corona and the Titan hydrogen torus model should provide a 
most interesting study. Preliminary modeling activities regarding the new 
data were initiated near the end of this project. The study of these data 
will be undertaken at AER in a new NASA funded continuation project for 
this research. 
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III. COMETARY ATMOSPHERES 


l . 


1 . Overview 

The three-year plan for this project in the area of cometary atmo- 
spheres is summarized in Table 4. The project has been very successful in 
meeting these objectives and has made significant advances in under- 
standing the physics and interrelationships in the structure and dynamics 
of the inner and extended comae of comets. 

In the first year, models for the hydrogen, carbon and oxygen comae 
were developed. In addition, model calculations for hydrogen and oxygen 
comae of Comet P/Encke, appropriate to the Pioneer Venus Orbiter 
Ultraviolet Spectrometer (PVOUVS) measurements acquired on 14 and 15 April 
1984 by A.I.F. Stewart, and also preliminary model calculations for the 
carbon coma observations of Comet Kohoutek acquired much earlier by Opal 
and Carruthers (1977) were performed (see the 1985 Annual Report) . As an 
outgrowth of the Comet P/Encke observations, a collaborative program was 
established with A.I.F. Stewart to analyze hydrogen data for Comet 
P/Giacobini-Zinner and H, 0, C and OH data for Comet P/Halley to be ac- 
quired by the PVOUVS in the second project year. 

In light of the significant progress accomplished in the first proj- 
ect year (see Table 4), efforts in the second project year were focused in 
three main areas: (1) analysis of PVOUVS hydrogen coma data for Comet 

P/Giacobini-Zinner where the effects of a time-variable solar wind and UV 
flux were also included, (2) performance of preliminary supporting model 
calculations and initial modeling analysis for PVOUVS observations of 
Comet P/Halley, and (3) generalization of the hydrogen comae model to in- 
clude collisions of non- thermally produced H atoms with the slower and 
heavier outflow species (i^O, OH, 0) and, in addition, testing the model 
with Comet Kohoutek observational data as preparation for analysis of the 
PVOUVS data for Comet P/Halley. Our analysis of the hydrogen coma data 
for Comet P/Giacobini-Zinner was published in Geophysical Research Letters 
(Combi, Stewart and Smyth 1986). The preliminary supporting model cal- 
culations for observation of Comet P/Halley are discussed in the 1986 
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COMETARY ATMOSPHERES: THREE-YEAR PLAN FOR MODELING ANALYSIS 
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Annual Report. Generalization of the hydrogen comae model to include col- 
lisions and the role of collisional thermal izat ion are discussed below. 

In the third project year, the two main areas of endeavor have been 
(1) assessment and preliminary analysis of the PVOUVS data for Comet 
P/Halley, and (2) completion of the hydrogen model including collisions 
and its successful employment in understanding the role of collisional 
thermal izat ion in the spatial morphology of the Lyman-a coma of comets. 
These two main areas are discussed below. 

2. Data and Preliminary Analysis for Comet P/Halley 

Coma data for Comet P/Halley were obtained by the 1985/1986 observa- 
tional program conducted by Dr. A.I.F. Stewart during the second project 
year using the ultraviolet spectrometer of the Pioneer Venus Orbiter. This 
program acquired a rich set of observations that are summarized in 
Table 5. As looted there, early pre-perihelion observations were made be- 
tween 28 December 1985 and 7 January 1986. Additional observations began 
again on 30 January (10 days pre-perihelion) and continued until 7 March 
(26 days post-perihelion). There is a large and consistent data set for 
coraetary 0, C and OH over this time period. For atomic hydrogen, most 
of the data provide radial brightness scans across the coma, but some of 
the H data (see Table 5) were taken so as to construct Lyman-a images. 

For the 0, C and OH comae, the data provide radial brightness scans which 
contain useful signal to noise information, but only within a few pixels 
of the nucleus. The 0, C and OH data are therefore valuable for deter- 
mining the relative abundance of these species near the nucleus. Modeling 
studies of the spatial structure of the coma can therefore be undertaken 
only for hydrogen. 

Reduction of the H coma data by Dr. Stewart has been slower than ini- 
tially anticipated because of two main factors: (1) uncertainties in 
the spacecraft pointing, and (2) uncertainties in the removal of the in- 
terplanetary Lyman-a background emission which is very spatially non-uni- 
form. A preliminary discussion of the H, 0 and C data has been given by 
Stewart (1987). 
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TABLE 5 

COMET P/HALLEY OBSERVATIONS FROM THE PIONEER VENUS ORBITER 


1985 


1986 


Date of 

Total Daily 
Observing 


Observing Time oer Species 

Chr} 






Observation 

Time (hr! 

H 

O 

C 

OH 

OTHER 

28 December 

15 

8 

7 


. 


29 December 

15 

7 

8 

_ 

- 

_ 

30 December 

16 

6 

1 

9 

- 

- 

31 December 

16 

5 

3 

- 

8 

- 

1 January 

13 

9 

- 

4 

- 

- 

2 January 

20 

4 

15 

1 

- 

- 

3 January 

16 

4 

4 

8 

- 

- 

4 January 

15 

4 

3 

- 

8 

- 

5 January 

14 

4 

10 

- 

- 

- 

6 January 

14 

7 

- 

7 

- 

- 

7 January 

7 

3 

4 

- 

- 

- 

30 January 

3 

3 

- 

- 

- . 

- 

31 January 

16 

10 

2 

2 

2 

- 

1 February 

11 

5 

2 

2 

2 

- 

2 February* 

18 

18 

- 

- 

- 

- 

3 February* 

19 

19 

- 

- 

- 

- 

4 February* 

19 

13 

1 

2 

1 

2 

5 February* 

19 

17 

- 

- 

- 

2 

6 February* 

19 

17 

- 

- 

- 

2 

7 February 

15 

9 

1 

1 

1 

3 

8 February 

11 

7 

2 

2 

- 

- 

9 February! 

16 

8 

2 

2 

2 

2 

10 February 

11 

8 

1 

2 

- 

- 

11 February 

16 

9 

1 

2 

1 

3 

12 February 

16 

9 

1 

6 

- 

- 

13 February 

19 

12 

6 

1 

- 

- 

14 February 

15 

10 

5 

- 

- 

- 

15 February* 

20 

20 

- 

- 

- 

- 

16 February 

16 

10 

4 

2 

- 

- 

17 February 

16 

9 

3 

2 

2 

- 

18 February 

18 

11 

2 

2 

3 

- 

19 February 

12 

9 

1 

2 

- 

- 

20 February 

19 

12 

2 

2 

2 

1 

21 February 

20 

13 

2 

2 

2 

1 

22 February 

16 

13 

1 

1 

1 

- 

23 February 

19 

13 

2 

2 

2 

- 

24 February 

19 

13 

2 

2 

2 

- 

25 February 

21 

14 

3 

2 

2 

- 

26 February 

14 

6 

3 

3 

2 

- 

27 February 

19 

11 

2 

3 

3 

- 

28 February 

22 

14 

3 

2 

3 

- 

1 March 

18 

12 

2 

2 

2 

- 

2 March 

19 

11 

3 

2 

3 

_ 

3 March 

21 

13 

2 

3 

3 

- 

4 March 

21 

13 

3 

2 

3 

- 

5 March 

15 

8 

2 

2 

3 

- 

6 March 

20 

12 

3 

2 

3 


7 March 

2 

- 

1 

1 

- 

- 


*Image data acquired on these days for hydrogen 
{Perihelion 
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PVOUVS Lyman-a data for Comet P/Halley for 7 days in the December 
1985 - March 1986 period have been received to date from Dr. Stewart for 
preliminary analysis. The background Lyman-a emission from these scan 
data has been only approximately removed. The scan data on each date that 
includes the nucleus have been selected for modeling analysis. Seven 
scans were thus constructed and fit by the Particle Trajectory Model 
(PTM) . The scan on 9 February 1986 and its model fit are shown in 
Figure 2. The observing dates for the seven scans and their H 2 O produc- 
tion rate as determined by the model fits are summarized in Table 6. A 
complete analysis of the Lyman-a scan and image PVOUVS data for hydrogen 
will be undertaken in a new project that will continue this research work. 

3. The Role of Collisional Thermal izat ion in the Coma 

The particle trajectory model (PTM) for cometary hydrogen was gene- 
ralized in the second year to include elastic collisions between non- 
thermal (-8 km sec'^, -20 km sec"^, and -30 km sec) hydrogen atoms (pro- 
duced by photodissociation of H 2 O and OH) and the heavy molecules and 
atoms (H 2 O, OH, 0) that flow radially outward from the nucleus. For 
comets with large production rates and small perihelion distances such as 
Comet Kohoutek and possibly Comet P/Halley, these collisional processes in 
the inner coma have been shown to be important in determining the velocity 
dispersion and hence spatial distribution of H atoms observed in Lyman-a 
emission of the outer coma. 

A paper entitled "The Role of Hydrogen Thermalization in Shaping the 
Lyman-a Coma", which briefly discussed the physics of the inner coma and 
its effect on the outer H coma, was presented (in poster format) at the 
"Symposium on the Diversity and Similarity of Comets" (6-9 April 1987, 
Brussels , Belgium) . This paper was also published by ESA (Combi and Smyth 
1987). 

The PTM used in the above paper is but one of a more general class of 
models for cometary gases. Addressing this larger subject, a major model 
documentation paper entitled "Monte Carlo Particle Trajectory Models for 
Neutral Cometary Gases. I. Models and Equations" and a companion paper en- 
titled "Monte Carlo Particle Trajectory Models for Neutral Gases. II. 
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INTENSITY (KILORAYLEIGHS) 


1000.00 


100.00 


DAY 40.440 


DISTANCE ALONG THE SCAN <10**6 KM) 


Figure 2 

Model Analysis of Pioneer Venus Lyman-a Observations of Comet P/Halley 

The points are a nucleus-crossing swath across the coma of Comet 
Halley recorded by the Pioneer Venus Orbiter UV spectrometer on 
9 February 1986. The line shows the PTM analysis. The water 
production rate is 1.74 x 10^® s ” * on that day. 







Table 6 


WATER PRODUCTION RATE FOR COMET P/ HALLEY FROM PIONEER VENUS ORBITER 


Date of 
Observation 

Heliocentric 
Distance (AU) 

H 2 O Production Rate 
(1030 molecules sec“ 

1985 30 December 

1.022 

0.41 

1986 1 February 

0.614 

1.37 

4 February 

0.598 

1.40 

9 February 

0.587 

1.74 

13 February 

0.594 

1.65 

18 February 

0.620 

1.61 

7 March 

0.802 

1.02 
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Spatial Morphology of the Lyman-a Coma" have been completed and published 
(Combi and Smyth, 1988 a,b). Both papers are included in the Appendix of 
this report. The first paper documents this new class of Monte Carlo par- 
ticle trajectory models which are based upon physical processes of atoms 
and molecules in the coma. This model provides a new level of realism in 
predicting the spatial distribution of observed species in cometary comae 
and in exploring the basic physics of the transition zone between true 
fluid- flow and free molecular flow. The second paper successfully applies 
the general model of the former paper to simulate the hydrogen Lyman-a 
image for Comet Kohoutek and includes near perihelion one image repre- 
senting an extreme case for collisional thermalization. The success of 
the model in this application is particularly noteworthy. 
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